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Abstract.  Published  differential  cross  section  data  for  heavy  particle  collisions  between  xenon  ions  and  neutral  xenon  has 
been  incorporated  into  plasma  simulations  for  electric  propulsion  modeling.  A  fit  has  been  made  to  the  published  data  in 
order  to  estimate  the  relative  contribution  from  charge  exchange  and  elastic  collisions  and  to  reduce  the  computational 
cost  of  utilizing  the  differential  cross  section  in  existing  numerical  models.  Since  the  published  profiles  do  not  include 
scattering  data  near  0  degrees,  the  differential  cross  section  was  assumed  to  be  constant  at  low  angles.  The  angle  at  which 
the  differential  cross  section  was  assumed  to  transition  from  the  constant  profile  to  the  fit  was  chosen  such  that  the 
differential  cross  section  integrated  to  the  published  total  cross  section  value  for  xenon  scattering.  In  order  to  make  the 
resulting  differential  scattering  curve  generally  applicable  to  other  types  of  collisions  with  dissimilar  collision  partners, 
the  profile  was  converted  from  the  laboratory  frame  into  center  of  mass  coordinates.  Each  time  a  scattering  event  was 
determined  to  take  place  in  the  electric  propulsion  modeling  codes,  a  scattering  angle  of  the  incident  particle  was  chosen 
using  a  cumulative  distribution  function.  The  behavior  of  the  target  particle  was  determined  using  conservation  of  energy 
and  momentum. 

Keywords:  Plasma  Simulation,  Monte  Carlo  Methods,  Electric  Propulsion. 
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INTRODUCTION 

Large  angle  scattering  of  primary  ions  emitted  from  electric  propulsion  devices  is  a  principal  area  of  concern  for 
spacecraft  integration.  Charge  exchange  or  elastic  collisions  between  high  velocity  ions  and  the  slow  plume  of 
neutral  gas  trailing  the  spacecraft  may  result  in  contamination  through  sputtering  and  deposition.  Until  recently,  the 
energy  and  angular  dependence  of  the  products  of  these  types  of  collisions  were  unknown  for  the  xenon  propellant 
commonly  used  in  electric  propulsion  devices  such  as  Hall  thrusters.  Therefore,  in  the  past,  simple  scattering  models 
have  been  employed  such  as  charge  exchange  collisions  free  of  momentum  transfer  and  elastic  collisions  resulting  in 
isotropic  scattering. 

Recent  measurements  taken  at  Hanscom  Air  Force  Base  using  a  guided-ion  beam  experiment  have  resulted  in  the 
calculation  of  ion-atom  interaction  potentials  for  both  single  and  doubly  charged  xenon  ions  [1],  These  interactions 
potentials  were  used  to  compute  the  absolute  differential  cross  section  at  a  typical  Hall  thruster  primary  ion  energy 
of  270  eV  per  unit  charge.  The  new  differential  cross  section  data  has  been  combined  with  previous  measurements 
of  total  charge  exchange  cross  sections  and  incorporated  into  the  Draco  module  of  the  electric  propulsion  plume 
model.  Coliseum,  using  a  Monte  Carlo  collision  technique  [2],  For  computational  speed,  a  simple  log-log  form  is  fit 
to  the  differential  cross  section  data  and  a  cumulative  distribution  function  is  used  to  randomly  choose  the  post- 
collision  scattering  angles.  Post-collision  velocities  are  determined  based  on  conservation  of  energy  and  momentum. 
For  greater  flexibility  in  future  simulations  involving  unlike  collision  partners,  all  collision-based  computations  were 
done  in  the  center  of  mass  reference  frame  and  converted  to  a  laboratory  frame.  This  work  discusses  the  details  of 
implementation  for  the  new  collision  model  in  the  existing  electric  propulsion  code. 
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DIFFERENTIAL  CROSS  SECTION 


In  Reference  1,  the  authors  have  experimentally  measured  the  differential  cross  section  of  Xe+  +  Xe  and  Xe2+  + 
Xe  at  laboratory  ion  energies  between  5  and  40  eV.  They  have  used  their  results  to  derive  interaction  potentials 
which  they  have  then  used  to  calculate  the  absolute  differential  cross  section  at  energies  relevant  to  the  electric 
propulsion  device  known  as  the  Hall  thruster.  The  results  presented  here  are  based  on  their  calculation  of  the 
differential  cross  section  at  an  ion  energy  per  unit  charge  of  300  V.  In  order  to  make  this  data  more  computationally 
tractable  in  existing  Monte  Carlo  collision  methods,  the  following  three  steps  have  been  performed: 

1.  Fitting  a  curve  to  published  differential  cross  section  data. 

2.  Determining  a  cutoff  angle  below  which  to  assume  differential  cross  section  is  constant. 

3.  Converting  from  the  laboratory  reference  frame  to  the  center  of  mass  reference  frame. 

In  order  to  reduce  the  computational  cost  of  choosing  a  scattering  angle,  a  fit  has  been  used  to  characterize  the 
differential  cross  section  for  each  ion  charge  state.  Noting  that  the  published  data  is  approximately  linear  on  a  log- 
log  scale  the  following  form  has  been  selected: 


d<J 

dQ. 


=  0A-ellOB-el  +(9O-0)A~cllO 


A  ct  -t  f\B  ct 


LAB 


(i) 


where  6?  is  the  laboratory  scattering  angle  in  degrees,  da/dfl  is  the  differential  cross  section,  and  the  subscript ,/./)/!  ’ 
denotes  the  laboratory  reference  frame.  The  above  form  has  the  added  advantage  that  it  separates  the  contributions 
of  both  elastic  collisions  and  charge  transfer  collisions  to  the  total  differential  cross  section.  The  optimized 
coefficients  are  listed  in  Table  1  and  a  comparison  between  the  published  data  and  the  curve  fit  is  shown  in  Figure  1 
for  each  charge  state. 
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Angle  [Degrees] 
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FIGURE  1.  Comparison  between  published  differential  cross  section  data  and  optimized  curve  fit  for  (a)  Xe+  +  Xe 

and  (b)  Xe2+  +  Xe. 


TABLE  1.  Coefficients  used  in  curve  fit  to  published  data. 


Coefficients 

Xe+  +  Xe 

Xe2+  +  Xe 

A 

el 

-2.02 

-2.34 

B 

el 

3.24 

3.17 

6 

el 

0.0045 

0.059 

A 

ct 

-1.098 

-1.43 

B 

ct 

1.53 

1.47 

6 

ct 

0.0022 

0.143 

The  differential  cross  section  data  in  Ref.  [1]  is  only  provided  between  the  angles  of  0.3  and  85.95  degrees.  Since 
a  significant  amount  of  scattering  occurs  between  0  and  0.3  degrees  and  between  89.95  and  90  degrees,  the 
differential  cross  section  in  this  region  must  be  specified.  Note  that  the  form  of  the  curve  fit  specified  by  Eqn.  (1) 
diverges  at  0  and  90  degrees.  Therefore,  the  differential  cross  section  is  assumed  to  reach  a  constant  maximum  value 
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at  extreme  angles.  The  angle  at  which  the  form  transitions  from  Eqn.  (1)  to  a  constant  value  is  given  by  0_el  for  the 
elastic  contribution  and  0_ct  for  the  charge  transfer  contribution.  These  angles  are  separately  chosen  such  that  the 
total  cross  section  is  given  by  the  published  charge  transfer  cross  section  specified  in  Ref.  [3],  namely  53.6  A2  for 
Xe+  +  Xe  and  23.7  A2  for  Xe2 1  +  Xe.  The  total  cross  section,  alol,  is  calculated  using  the  following  formula: 


a„ 


=  1 


da 
0  dQ. 


2  7i  s'm(d raci)dOl 


rad 


(2) 


LAB 


where  6rad  is  the  laboratory  scattering  angle  in  radians.  When  integrating  over  the  charge  transfer  cross  section,  the 
angle  in  the  integral  is  replaced  with  (n/2  -  drad)  since  the  laboratory  scattering  angle  in  this  case  the  recoil  angle  of 
the  target  particle  which  was  initially  neutral.  The  resulting  values  of  0_el  and  6_ct  are  given  in  Table  1.  Figure  2 
demonstrates  the  constant  assumed  differential  cross  section  for  the  elastic  and  charge  transfer  scattering  of  Xe+  + 
Xe.  Note  that  the  elastic  total  cross  section  has  been  assumed  to  be  equal  to  the  charge  transfer  total  cross  section  for 
each  ion  charge  state. 


(a)  (b) 

FIGURE  2.  Assumed  form  of  Xe+  +  Xe  differential  cross  section  for  (a)  elastic  and  (b)  charge  transfer  scattering. 


For  ease  of  conserving  momentum  and  energy,  as  well  as  to  generalize  the  resulting  profiles  to  unequal  collision 
partners,  the  differential  cross  section  was  converted  from  laboratory  to  center  of  mass  coordinates  using  the 
following  formulas  for  collision  partners  of  equal  weight: 


da 

da 

1 

dD. 

cm  dCl 

LAB  4C0S# 

(3) 


%  =  20  (4) 

where  %  is  the  center  of  mass  scattering  angle  in  degrees  and  the  subscript  „CV/’  indicates  center  of  mass.  Note  that 
for  charge  transfer  collisions,  9  is  replaced  with  90  -  0. 


In  terms  of  the  center  of  mass  scattering  angle,  %,  rather  than  the  solid  angle,  Q,  the  differential  cross  section  can 
be  expressed  as: 

da  da 

- = -  2^-siny  (5) 

dx  dnCM 

Therefore,  using  the  trigonometric  identity,  sitl(2f?)  =  2  COS(t9)  S  i  tl  ( (9) ,  the  differential  cross  section  for  elastic 
collisions  is  given  by: 
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(6) 


(7) 


Several  different  methods  are  employed  to  incorporate  heavy-particle  collisions  into  the  existing  electric 
propulsion  models.  The  first  step  in  each  method  is  determining  which  particles  in  the  simulation  should  undergo  a 
collision  at  each  time  step. 

In  the  existing  electric  propulsion  codes,  the  decision  to  perform  a  collision  is  determined  by  testing  every  ion 
particle  individually.  First,  a  collision  partner  is  chosen  randomly  either  by  constructing  a  random  particle  using 
fluid  background  properties  or  by  choosing  a  particle  randomly  from  within  the  same  cell.  Using  the  properties  of 
the  ion  being  considered  and  the  randomly  chosen  neutral  partner,  the  probability  of  a  collision  is  calculating  using 
the  following  formula: 

P  =  1  -  exp  (-  n  o\cr  |  At )  (8) 


where  P  is  the  probability  of  collision,  n  is  the  number  density  of  the  target  species  (neutrals  in  this  case),  c,  is  the 
relative  velocity,  and  At  is  the  timestep.  This  probability  is  compared  to  a  random  number  to  determine  whether  or 
not  a  collision  takes  place. 

Once  it  has  been  determined  that  a  collision  is  statistically  likely  to  occur  and  the  collision  partners  have  been 
chosen,  the  following  steps  are  performed  to  execute  the  collision. 

1.  Choose  center  of  mass  scattering  angles  using  statistical  sampling. 

2.  Assuming  pre-collision  relative  velocity  is  aligned  along  a  principle  direction,  calculate  post-collision 
relative  velocity. 

3.  Rotate  coordinate  axes  to  account  for  actual  direction  of  pre-collision  relative  velocity. 

4.  Convert  from  center  of  mass  to  laboratory  reference  frame. 

In  carrying  out  a  collision,  the  first  step  is  to  use  statistical  sampling  to  determine  the  polar  and  azimuthal 
scattering  angles.  For  convenience,  it  is  easiest  to  align  the  pre -collision  relative  velocity  along  a  principle  direction, 
such  as  the  z-axis  as  shown  in  Fig.  3  (a).  Due  to  axisymmetry,  the  azimuthal  scattering  angle  is  chosen  randomly 
from  0  to  2 n.  The  polar  scattering  angle  is  chosen  using  the  differential  cross  sections  specified  by  Eqns.  6  and  7 
depending  on  the  collision  type.  For  each  type  of  collision,  a  comparison  function  is  used  to  select  a  random  angle 
from  the  distribution.  Based  on  the  form  of  the  differential  cross  section,  the  comparison  function  chosen  was: 


f'iz) 


\eA  for  x<e 

\xA  f°r  x>° 


(9) 


where  9  and  A  are  the  coefficients  specific  to  each  type  of  collision  as  listed  in  Table  1.  In  order  to  select  a  number 
randomly  from  this  comparison  function,  a  cumulative  distribution  function  is  used.  The  polar  scattering  angle,  is 
chosen  such  that: 
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(10) 


\fl{a)da 
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J /,(«}/« 


=  £ 


0 

where  is  a  random  number  between  0  and  1.  An  acceptance-rejection  technique  is  then  used  to  determine  whether 
or  not  to  select  this  value  of  %  or  pick  a  new  value.  Based  on  the  form  of  the  differential  cross  section,  the  value 
selected  is  retained  with  probability  sin(^ /  2 )  for  elastic  collisions  and  COS (jj^/2)  for  charge  transfer  collisions. 


Once  azimuthal  scattering  angle,  (|>,  and  polar  scattering  angle,  %,  have  been  determined,  the  post-collision 
relative  velocity  in  the  z-aligned  reference  frame  is  computed  using  the  following  formulas: 

cfx  =  c;  sin  x  cos  <j> 

czt*  -  czr  sin x sin (f>  (11) 

cz*  =c;  cos  x 


where  the  subscript  ,/•’  denotes  relative  velocity,  the  superscript  denotes  a  z-axis  aligned  reference  frame,  and  the 
superscript  denotes  a  post-collision  velocity. 


Since  the  above  formulas  assume  that  the  initial  relative  velocity  was  aligned  along  a  principle  direction  (the  z- 
axis),  a  rotation  must  be  performed  on  the  resulting  relative  velocity  vector  to  account  for  the  actual  orientation  of 
the  initial  relative  velocity  vector.  Assuming  that  the  polar  and  azimuthal  rotation  angles  are  given  by  P  and  y 
respectively,  as  shown  in  Fig.  3  (b),  the  following  transformation  yields  the  final  relative  velocity  vector: 


—  * 


p = tan~’  Ucl  +c 


/o 


cos  y 

-sin  y 

0" 

cos  1 3 

0 

sin  P 

sin  y 

cosy 

0 

0 

1 

0 

0 

0 

1 

-sin  P 

0 

COS  P 

(12) 

(13) 


(14) 


Once  the  post-collision  relative  velocity  has  been  determined  in  the  center  of  mass  frame,  the  final  velocity 
vectors  of  each  particle  in  the  laboratory  reference  frame  are  calculated.  If  the  collision  type  is  elastic,  these 
velocities  are  directly  applied  to  the  source  and  target  particles.  In  the  case  of  charge  exchange  collisions,  the  final 
velocity  of  the  source  particle  is  applied  to  the  target  particle  and  the  final  velocity  of  the  target  particle  is  applied  to 
the  source  particle  to  reflect  the  transfer  of  charge  that  accompanies  this  type  of  collision. 


TEST  CASE 


In  order  to  test  the  functionality  of  the  collision  routine  implemented  in  the  electric  propulsion  plume  model, 
COLISEUM,  a  simple  test  case  was  constructed  in  order  to  verify  that  the  imposed  differential  cross  section  is 
correctly  utilized  by  the  code.  The  test  case  geometry,  illustrated  in  Fig.  3(a),  consists  of  a  cylindrical  beam 
impinging  on  a  rectangular  target  box  of  neutral  gas.  The  total  domain  size  of  the  simulation  is  10  cm  x  10  cm  x  2 
cm  with  a  cell  spacing  of  1  mm,  resulting  in  200,000  cells  in  total.  The  injector  surface  is  a  10  mm  diameter  circle 
with  a  flow  rate  of  Xe'  of  2.78e-8  kg/s.  The  directed  energy  of  the  beam  is  300  eV,  corresponding  to  a  velocity  of 
21144  m/s,  with  a  thermal  energy  component  of  100K.  Under  these  conditions,  the  density  of  the  beam  is 
approximately  8el6  m'3.  The  target  is  12  mm  x  12  mm  x  2  mm  in  size  with  an  initial  density  of  2E19  m"3 
represented  by  576,000  macroparticles  of  weight  1E7  neutrals  per  macroparticle.  Based  on  these  conditions,  the 
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collisional  mean  free  path  is  93.45  mm  for  singly  charged  ions  and  212  mm  for  doubly  charged  ions.  Since  the 
dimension  of  the  target  box  parallel  to  the  beam  direction  is  2  mm,  it  is  highly  unlikely  that  a  source  particle  would 
undergo  more  than  one  collision. 


FIGURE  3.  Coliseum  test  case  (a)  geometry  and  (b)  comparison  of  calculated  and  imposed  differential  cross 

section. 

The  position  and  velocity  information  of  all  particles  crossing  a  sample  box  of  size  16  mm  x  16  mm  x  6  mm 
around  the  target  is  stored  for  post-processing.  The  velocity  information  of  the  sampled  ions  is  used  to  calculate  a 
scattering  angle  for  each  particle  which  is  binned  to  generate  the  calculated  differential  cross  section.  This  profde  is 
normalized  to  the  published  differential  cross  section  at  45  degrees.  As  shown  in  Fig.  3(b),  the  generated  differential 
behavior  is  consistent  with  the  imposed  differential  cross  section. 

SUMMARY 

Using  published  differential  cross  section  data  for  xenon  ions  colliding  with  neutral  xenon,  a  technique  has  been 
developed  to  utilize  this  data  in  an  existing  electric  propulsion  code,  COLISEUM.  The  data  has  been  pre-processed 
to  be  consistent  with  published  total  cross  section  values  and  to  isolate  the  contributions  due  to  elastic  and  charge 
exchange  collisions,  as  well  as  for  computational  speed. 

ACKNOWLEDGMENTS 

The  authors  wish  to  thank  Y u-Hui  Chiu  and  Rainer  Dressier  for  providing  the  raw  differential  cross  section  data 
at  300  eV  and  for  their  continued  assistance  in  applying  their  measurements  to  electric  propulsion  modeling. 

REFERENCES 

1.  Y.-H.  Chiu,  R.A.  Dressier,  D.J.  Levandier,  C.  Houchins,  and  C.Y.  Ng,  “Large-angle  xenon  ion  scattering  in  Xe-propelled 
electrostatic  thrusters:  differential  cross  sections,”  Journal  of  Physics  D:  Applied  Physics,  41,  165503,  2008. 

2.  L.  Brieda,  R.  Kafafy,  J.  Pierru,  and  J.  Wang,  “Development  of  the  DRACO  Code  for  Modeling  Electric  Propulsion  Plume 
Interactions,”  Proceedings  of  the  40th  Joint  Propulsion  Conference,  AIAA-2004-3633,  2004. 

3.  J.  S.  Miller,  S.H.  Pullins,  D.J.  Levandier,  Y.-H.  Chiu,  and  R.A.  Dressier,  “Xenon  charge  exchange  cross  sections  for 
electrostatic  thruster  models,”  Journal  of  Applied  Physics,  91,  3,  2002. 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


THE  AIR 


FORCE  R 
LEAD  I 


ESEARCH  LABORATORY 
DISCOVER  I  DEVELOP  I  DELIVER  ir 


DSMC  Implementation  of  Experimentally-Based 
Xe+  +  Xe  Differential  Cross  Sections  for  Electric 

Propulsion  Modeling 

11-15  July  10 


Michelle  Scharfe,  Justin  Koo, 

Gregory  Azarnia 

EP  Modeling  and  Simulation 
Propulsion  Directorate 
Air  Force  Research  Laboratory 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Overview 


*  Introduction 

*  Differential  Cross  Section 

*  Monte  Carlo  Methods 

*  Testing  Collision  Routines 

*  Summary 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


2 


Introduction 


9  Motivation  for  Electric  Propulsion 
Modeling 

—  Understand  device  physics 
—  Thruster  optimization 
—  Spacecraft  integration 

9  Collision  Models 


—  Collisions  important  for  understanding 
sputtering  and  deposition 

—  Before:  Isotropic  scattering  and  simple 
charge  transfer 

—  Now:  Based  on  published  differential  cross 
sections  [i] 

[1]  Chiu,  Dressier,  et  al.  “Large-angle  xenon  ion  scattering 
in  Xe-propelled  electrostatic  thrusters:  differential  cross 


sections,”  JOP  D,  41, 


1  65^3>u^QQ8Approved 


for  public  release;  distribution  unlimited 


3 


[js/y]  fJp/fp 


Differential  Cross  Section: 
Fitting  a  Curve  to  Published  Data 


*  Xe  ion/neutral  differential  cross  section  published  by  Chiu,  et  al. 

*  First  step:  Fit  a  curve  to  published  data  (300  eV) 

•  Increased  computational  speed 

•  Allows  separation  of  elastic  and  charge  transfer  contributions 


da 


dfl 


=  0A-ellO 


\B 


LAB 


d  +{9Q-0)A-a\0B-a 


Xe+  +  Xe 


Angle  [Degrees] 


Angle  [Degrees] 
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Differential  Cross  Section: 
Determining  Cutoff  Angle 


*  Published  data  does  not  contain  information  below  0.3  degrees 
or  above  89.95  degrees. 


Second  step:  Use  published  total  cross  section  values  to 
determine  cutoff  angle[2]:  da 


I 


dd 


2*- sin  (dmd)d0rad 


LAB 


Anqle  [Deqreesl 


Xe+  +  Xe  Charge  Transfer  Scatterinq 


[2]  Miller,  Pullins,  et  al.,  “Xenon  charge  exchange  cross  sections  for  electrostatic  thruster 
models,”  Journal  of  Applied  Physics,  91,  3,  2002. 
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Differential  Cross  Section: 
Converting  to  Center  of  Mass  Frame 


•  Easier  to  perform  rotations  in  center  of  mass  frame 

9  Also,  makes  differential  cross  section  more  meaningful  when 
applied  to  dissimilar  collision  partners 

•  Third  step:  Convert  differential  cross  section  to  center  of  mass 
coordinates 


da  da  1 


d£l  CM  d£l  LAB  4  cos  6 


Z  =  26> 


where  0  is  the  lab  frame  scattering  angle  and  x  is  the 
center  of  mass  frame  scattering  angle. 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


6 


Differential  Cross  Section: 
Final  Form 


•  In  order  to  use  differential  cross  section  to  choose  scattering 
angle,  convert  solid  angle  to  center  of  mass  scattering  angle. 

da  dci 

- = -  2/rsin  x 

dx  dCl  CM 


Final  form  of  differential  cross  section: 

A  el 


dcr 

dx 


X  I  e i 


;rsm|  — 


{6  _el)A-elVdB~el  7ism{6  _el) 


for  %  >  2d  _el 
for  x  <  20 _el 


da 

dx 


10* 


ct 

71  COS 


[0  _  ct  )A~ct  1 0B-Ct  n  cos(<9  _ct) 


for  X-2 0 _ct 
for  x  <20  _ct 
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Monte  Carlo  Methods: 
Choosing  Collision  Partners 


*  Several  mechanisms  for  choosing  collision  partners 

•  Calculate  the  probability  of  each  source  particle  colliding  based  on 


target  density  (COLISEUM) 


•  Use  a  no-time-counter  method  to  randomly  pick  appropriate 
number  of  pairs  (HPHall). 

*  Once  collision  partners  have  been  chosen: 

•  Azimuthal  scattering  angle  is  chosen  uniformly:  0  <  4>  <  2n 

•  Use  differential  cross  section  equations  to  choose  polar  scattering 
angle. 

•  Apply  scattering  angles  to  initial  velocity  vectors  to  determine  post¬ 
collision  velocities. 
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Monte  Carlo  Methods: 
Choosing  Scattering  Angles 


*  To  apply  statistical  sampling  from  differential  cross  section,  first 
choose  a  comparison  function: 


9A  for  X-9 
\xA  far  X>  9 


•  Then  randomly  pick  an  angle,  %,  from  comparison  function  using 
cumulative  distribution  function  method. 


Where  i;  is  a  random  number:  0<«^<1 


*  Use  an  acceptance-rejection  technique  to  decide  whether  or  not 
to  retain  the  selected  angle  or  choose  again. 
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Monte  Carlo  Methods: 
Calculating  Post-Collision  Relative  Velocity 


*  For  simplicity,  assume  that  the 
pre-collision  relative  velocity  was 
aligned  with  a  principal-axis, 
such  as  z. 

*  Calculate  the  post-collision 
relative  velocity  in  this  rotated 
reference  frame: 

cll  =  cr  sin  X  COS  <]) 
czr*  =  czr  sin  x sin  <f> 
cz*  =  cl  cos  X 
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Monte  Carlo  Methods: 
Rotating  Coordinate  Frame 


•  To  account  for  altered  reference  frame,  first  calculate  rotation 
angles  to  align  pre-collision  relative  velocity  vector  with 
principle  direction: 


P  =  tair‘  W  Cl+Cy  ) 

fC  \ 

y  =  tan  1  — 

\CX  J 

•  Rotate  post-collision  relative  velocity 
vector  according  to  rotation  angles: 


cos  y 

-sin ;/ 

0" 

cos  p 

0 

sin  p 

sin^ 

cos  y 

0 

0 

1 

0 

0 

0 

1 

-sin  p 

0 

COS  P 
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Monte  Carlo  Methods: 
Converting  to  Laboratory  Frame 


Lastly,  convert  back  to  laboratory  reference  frame  from  center  of 
mass  reference  frame: 


C 


m 


—  c  + 

source  m 


target 


C, 


TYL  +  TYL 

L  source  ''target 


C 


m 


target 


source 


C, 


TYL  +  YYL 

source  v target 


•  In  the  case  of  charge  transfer  collisions,  the  source  and  target 
particle  identities  are  exchanged  after  the  new  velocities  have 
been  calculated. 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


12 


COLISEUM  Test  Case 


*  COLISEUM  is  an  electric  propulsion 
plume  model  using  hybrid  techniques 
with  PIC  heavy  particles 

9  Generate  simple  case  in  COLISEUM  to 
test  collision  routines 

•  Injector:  Xe+,  10  mm  diameter,  300  eV,  100 
K,  2.78e-8  kg/s,  8E16  m-3 


•  Target  cell:  12  mm  x  12  mm  x  2  mm, 
576,000  macroparticles,  2E19  rrr3 

•  Mean  Free  Path:  93.45  mm 

9  Sample  ions  crossing  a  box  (16  mm  x 
16  mm  x  6  mm)  around  target  cell 

9  Angular  distribution  of  sampled  ions  is 
consistent  with  imposed  differential 
cross  section 


12  mm 
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12  mm 


Summary 


*  Differential  cross  section  data  for  collisions  between  xenon  ions 
and  neutral  xenon  has  recently  been  published. 

*  This  new  data  allows  for  more  accurate  predictions  of  post¬ 
collision  products. 

*  The  published  differential  cross  section  has  been  simplified  and 
pre-processed  for  use  in  existing  electric  propulsion  models. 

*  Simple  test  case  in  electric  propulsion  plume  code,  COLISEUM, 
verifies  proper  implementation. 

9  Using  the  new  collision  data  allows  for  more  accurate  predictions 
of  the  interactions  between  the  high  velocity  ions  leaving  the 
device  and  the  slow-moving  neutral  plume  of  exhaust. 

•  Provides  more  accurate  calculation  of  high  angle  fluxes 

•  Improves  spacecraft  integration  modeling 
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